Two-terminal conductance fluctuations in the integer quantum Hall regime 
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Motivated by recent experiments on the conductance fluctuations in mesoscopic integer quantum 
Hall systems, we consider a model in which the Coulomb interactions are incorporated into the 
picture of edge-state transport through a single saddle-point. The occupancies of classical localised 
states in the two-dimensional electron system change due to the interactions between electrons when 
the gate voltage on top of the device is varied. The electrostatic potential between the localised 
states and the saddle-point causes fluctuations of the saddle-point potential and thus fluctuations of 
the transmission probability of edge states. This simple model is studied numerically and compared 
with the observation. 



PACS: 73.20, 73.40.Hm. 

I. INTRODUCTION 

Mesoscopic sample-to-samplc fluctuations in the con- 
ductance, G, are universal in the sense that the devi- 
ation from its average, (SG) 2 — (G — (G)) 2 , is alwa; 
about the order of (e 2 /h) 2 at very low temperaturetr 
In zero or weak external magnetic fields, this behaviour 
can be understood in termsr-iof the quantum phase co- 
herence of diffusive electronsa. This picture, however, is 
not expected to work in the quantum Hall regime since 
the trajectories of electrons are drastically deformed in 
strong magnetic fields. In this situation, the existence 
of extended states at sample edges makes the transport 
properties rather different from the weak-field caseo. In- 
deed, magnetoconductoace measurements on mesoscopic 
multiterminal devicesBEl seem to agree with this expec- 
tation. To be consistent with these observations, explicit 
inclusion of ttie-edge-channel effects in theories seem to 
be necessaryE2nl3. 

The appearance of the bulk delocalised states when 
the Fermi level is between two quantised Hall plateaus 
makes a mesoscopic system in strong magnetic field 
even more interestingtj. It is only recently, how- 
ever, that experimental data on the statistical be- 
haviour of the— .conductance in this regime has be- 
come availablella. Specifically, conductance fluctua- 
tions has been studied in mesoscopic n-channel Silicon 
Metal-OxidfcSemiconductor Field-Effect Transistors (Si 
MOSFETs)E£l with two terminals. In these systems, con- 
ductance is quantised (in units of e 2 /h) at some values 
of the gate voltage, V g , and the complication due to the 
Shubnikov-de Haas oscillation in multiterminal devices is 
absent in the two-terminal setting. It is found that the 
distribution of conductance below the first conductance 
plateau is almost uniform between G = and e /h. This 
is in clear contrast to the Gaussian distribution in the 
weak-field caseEfl. Theoretical studies have been carried 
out on models based on single-paiiicle interference effect, 
e.g. tight-binding—Hamiltonian 13, Chalker-Coddington 
network modelUlij, and so onqLpumerically and also 
a renormalisation group analysisEd. Some capture this 
feature of the conductance distribution. 



More recently, Cobden and collaborators demonstrate 
further that there is a great difference between conduc- 
tance fluctuations . in. the lowest few Landau levels and 
the weak field caseE3n3. Using the same sets of Si MOS- 
FET devices as described, they observed in strong fields 
strong correlations of fluctuation peaks and dips at differ- 
ent magnetic fields. This is achieved by studying explic- 
itly the evolution of fluctuation peaks and dips with mag- 
netic fields in different field regimes. Using a greyscale 
plot of dG I dV g as a function of B and V g , each bright or 
dark line can be viewed as the 'history' line of a peak or 
a dip of the conductance at different combinations of V g 
and B. The contrast between weak and strong magnetic 
fields is very clear (See Fig. 2 in Ref. 22). 

For weak magnetic fields, no pattern can be seen what- 
soever. At B >10T, we can see peaks and dips of conduc- 
tance fluctuations in a transition region between two suc- 
cessive conductance plateaus, which are wide and grey re- 
gions. The fluctuations 'move' in gate voltage, in straight 
lines as the magnetic field is varied. This pattern means 
that there exist strong correlations between fluctuation 
peaks and dips at different gate voltages at high magnetic 
fields. A closer look reveals that there are two types of 
lines, with different slopes, in a transition region. The 
slope of one set of lines is parallel to that associated with 
the centre of the plateau below, while the slope of the 
other set of lines is parallel to that associated with the 
centre of the plateau above. 

The strong contrast between the fluctuations on the 
Vg-B plots in weak and strong magnetic fields means that 
what occurs in two different field regimes is dramatically 
different. The straight lines appearing in the strong-field 
results makes it no longer appropriate to treat the_ahc=. 
nomenon as a single-particle one. Thus those modelsEjtiJ 
we mentioned before are obviously not compatible with 
strong correlations of conductance fluctuations at differ- 
ent fields. 

In this paper, we consider a model taking into ac- 
count the electron-electron interactions. The conduc- 
tance plateaus are understood in terms of the edge-state 
picture and the absence of backscattering between chan- 
nels at opposite edges of the sample. In the transition 
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region between plateaus, we assume that the paths fol- 
lowed by electron currents coming from two contacts (or 
terminals) percolate into the bulk as the chemical po- 
tential is increased and eventually become connected by 
tunneling through a single saddle point. In addition to 
the electrons in these edge states, we also take account 
of electrons in localised states. Conductance fluctuations 
arising from Coulomb interactions between both sets of 
electrons are treated in a simple way. 

As electrons are added to the system, the ones already 
in the localised states will rearrange due to the interac- 
tions. Therefore, the occupation of a given localised state 
can fluctuate several times between and 1 when the 
chemical potential, /i, is changed from — oo to +00. The 
total electrostatic potential between the localised states 
and the saddle point then also fluctuates as [i is varied. 
Because of the fluctuations of the saddle point potential, 
the tunneling of percolation paths fluctuates as a func- 
tion of the chemical potential. In our model, fluctuations 
of the conductance are purely due to these interaction 
effects. Numerical study of this model shows that the 
conductance between plateaus as a function of the chem- 
ical potential does indeed fluctuate significantly within 
the range between G = and e 2 /h. 

In the next section, we first develop our model, based 
on the transmission of edge-state channels through a sin- 
gle saddle point in the transition region. Fluctuations of 
a single potential saddle point are related to the Coulomb 
interaction between localised states in the bulk of the 
disordered 2DES. These localised states are treated us- 
ing the Efros-Shklovskii Coulomb glass model. Models 
at finite temperature T with localised states arranged on 
regular lattice points or random sites are then studied 
numerically by Monte Carlo methods and exact enumer- 
ation. Results for G versus V g at different disorder con- 
figurations and temperatures are obtained. After that, 
comparisons of our numerical results and experiments are 
discussed. We conclude with a summary. 



II. THEORETICAL MODEL 



MOSFETsME3€3, the bias voltage between two contacts 
is fixed, but the gate voltage of the metallic gate on top 
of the 2DES, therefore the Fermi energy, is varied. As 
the number of electrons is increased such that the Fermi 
level is between Landau levels E n < Ep < E n+ i, there 
are n edge channels at each boundary. In a high magnetic 
field with no inter-edge backscattering, the conductance 
of the system is thus quantised and given by G = ne 2 /h. 

In the transition region between two quantised 
plateaus, it is clear that the Fermi level is within a 
disorder-broadened Landau level. For clarity, let us con- 
sider what happens when the conductance is between the 
first and the second plateaus (see Fig. Q). At zero tem- 
perature, there are localised states hi-Jjhe bulk associated 
with closed local equipotential lincsE^, shown as dashed 
and dotted closed lines in Fig. ^. Apart from them, we 
can imagine that, in the bulk of the system, there are also 
directed extended states at the boundaries of the regions 
which are occupied by electrons and connected to one or 
the other contact. At some value of the chemical poten- 
tial, these extended states undergo large excursions into 
the bulk. As fi is increased further, the directed extended 
states approach each other more closely and tunneling 
between them occurs (shown as the dotted line in Fig. 
|l|) . Eventually, a second pair of edge states is formed in 
this way at the inner side near the system boundaries. 



FIG. 1. The formation of an edge-state channel pair as 
the gate voltage is increased in the transition regime, here 
between G = e 2 /h and 2e 2 /h. The lines with arrows represent 
the directed trajectories of electrons or holes in the magnetic 
fields. 



Before introducing the model, we note that there are 
different electron states in a disordered two dimensional 
electron system (2DES) in a high magnetic field. Edge 
state exists on an equi-energy line along the bound- 
ary of the sample connecting two contacts. The width 
of each edge-state channel is about the order of mag- 
netic length, Ib oc 1/y/B. Therefore, in a high mag- 
netic field the backscattering between edge states at op- 
posite boundaries can be neglected. Since each of the 
edge-state channels is now one-dimensional, we can ap- 
ply the two-terminal Landauer-Biittiker formula to cal- 
culate the conductance. The current injected into each of 
the edge states is proportional to the difference of chem- 
ical potentials between two contacts and also each edge- 
channel pair. In the context of the measurements on Si 



Let us consider the simplest tunneling process in which 
the two growing regions of electrons form a single poten- 
tial saddle point assumed to be located at about the cen- 
tre of the 2DES. The transmission probability 7y = |ty| 2 
for the incoming channel i and outgoing channel j is al- 
ready known, for the saddle-point potential of the form 
V(x, y) = V — (l/2)muo 2 x 2 + (l/2)mu) 2 y 2 with lu x and 04 
characterising the shape of the saddle-point potentials. 
It is given by 



TlJ SlJ 1 + exp(-7re„) 



(1) 



where e n = [E— E2(n+ 1/2) — V s -p]/Ex with E denoting 
the total energy of the electron and V s - V = Vq being the 
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bare potential strength. The two-terminal conductance 
is then given by G = (e 2 /h) £\ T u . 

In Eq. (0), E x and E 2 are in general complicated func- 
tions of lo c (cyclotron frequency) , lu x , and oj v . At high 
magnetic fields, when u c is much bigger than u) x and ui y , 
e n in Eq. ([!]) is simply a dimensionless measure of the 
energy of the electron's guiding-centre motion relative to 
Vq. For the edge state of the percolating region with a 
guiding centre energy Eq = Ep — 7iu; c (n+l/2), the trans- 
mission probability is one if Eq ^ Vq. In the other limit 
that E G <C Vo, the edge state is completely reflected. The 
two-terminal conductance as a function of Fermi energy 
at various magnetic-field-, strengths have been obtained 
numerically by Biittikerca. Interestingly, although these 
results are for a 2DES with a point-contact constriction 
or split-gates on top, the sape structure is observed in 
a 'macroscopic ' Si MOSFETEl. From this, it seems that 
considering only one potential saddle-point in the meso- 
scopic Si MOSFET where measurements were carried out 
may be reasonable. 

Within the framework of the edge-state picture de- 
scribed above, fluctuations of the transmission of perco- 
lating states lead to the conductance fluctuations. Here 
the transmission of electron currents depends solely on 
the potential energy at the saddle point. In general, this 
potential energy should have a contribution from the in- 
teractions. It is our purpose to relate the interactions 
to the conductance fluctuations. Electrons in our model 
are separated into categories: those in edge states and 
those in localised states. In principle, a full treatment-af 
two interactions between electrons is very complicatedc3. 
For simplicity, we neglect interactions amongst electrons 
in edge states, and just consider interactions between dif- 
ferent electrons in localised states and their effect on the 
saddle-point potential. Equally, for simplicity, wc neglect 
tunneling processes that mix the states we term localised 
with each other or with those we term edge states. This 
is correct in a finite sample in the semiclassical limit, 
when for a given value of the Fermi energy tunneling 
is important only at one saddle-point of the potential. 
More generally, one expects states at the plateau tran- 
sition to be delocalised by the tunneling processes we 
omit, at least within a single-particle description. Note 
that, although interactions between extended electrons 
are left completely untreated in the model, it turns out, 
as we shall present later, that this way of including the 
Coulomb interaction does indeed produce dramatic fluc- 
tuations of the saddle-point potential as the chemical po- 
tential of the system is varied. 

The model Hamiltonian describing a system of lo- 
calised states interacting with the Coulomb interaction 
is essentially the one which was studied first by Efros 
and ShklovskiiEZl. The same Hamiltonian produces a gap 
in the single-electron density of states at the Fermi level 
due to the Coulomb interaction. In the 'Coulomb-gap' 
system, electrons are strongly localised on a discrete set 
of sites due to the impurities. These quantum particles 
can thus be thought of as being in the regime where they 



behave classically. The Hamiltonian for localised states 
on N sites at positions denoted by i is 
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H = Yei(ni--) + -Y 



U 1^ e 2 (n 4 -I)K-I) 



N 1 
/i^(n 4 - -), 



(2) 



i=l 



where e< is the random site-energy, rij is the occupation 
number which can be either or 1, and r,j is the dis- 
tance between site i and j. The subtraction of 1/2 from 
Hi represents neutralising background charge, e and Sq 
are dielectric constants for the Silicon and the vacuum, 
respectively. The last term in the Hamiltonian is crucial 
to model the experimental variation in the chemical po- 
tential in the system as the gate voltage is varied. More 
explicitly, as the voltage of the gate electrode on top of 
the 2DES is varied, there will be electrons brought in or 
pulled out of the 2DES in the inversion layer in the Sil- 
icon. The occupancy in each of the localised states will 
then fluctuate as a function of /j,, due to the Coulomb 
interaction. Accordingly, the electrostatic potentials be- 
tween the N localised states and the saddle-point poten- 
tial will also fluctuate. Hence we have a fluctuating total 
saddle-point potential 
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with ri^ s -p representing the distance between site i and 
the saddle point, where h is the perpendicular distance 
between the inversion layer containing localised states 
and the centre of the saddle point. Again, the occupa- 
tion at site i is given by (rij — 1/2) to take into account 
background charges. Note that at finite temperature T, 
in Eq. (||) we need the thermal average of the occupation 
number, (rii), instead of rij at each site i. As a result, 
the transmission probability given by Eq. (Jj]) fluctuates 
as a function of the gate voltage. 



III. NUMERICAL RESULTS AND 
COMPARISONS WITH EXPERIMENTS 

In this section, we present the results obtained numer- 
ically at finite temperature. The 2DES is chosen to have 
a rectangular shape with length L and width W in units 
of lattice spacing. The localised states are chosen to be 
fixed at sites i arranged on a regular lattice with posi- 
tions given by (x,y). At fixed temperature T, we com- 
pute (n(x,y)), at each site i with coordinates (x,y) for 
each chemical potential /1. Note that, here and in all the 
following computation, we have set the electric charge, 
47T££o, and the Boltzmann constant ks to be one. For 
each configuration of impurities, a random energy is at- 
tached to each localised state. The random site-energy 
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e, at site i is some value in the interval [— W/2, W/2]. 
After (n(x,y)) is obtained, Eq. (||) determines the value 
of V s _ p at this [i and temperature T. The corresponding 
transmission probability is then given by Eq. ([!]), rewrit- 
ten here as T« = 1/1 + cxp[— (I4-p + c)/ i?o] , where c is a 
constant representing the energy .of the electron's guiding 
centre (or the equipotential line)c3; c, however, is always 
set to be equal to Vq in our numerical study. Eq here 
is an energy scale, chosen arbitrarily, which in principle 
can depend on the magnetic field and the characteristics 
of the saddle-point. The conductances at different chem- 
ical potentials are calculated in this way by employing 
the same set of random energies at N sites. For different 
sets of random site-energies which essentially represent 
different disorder realisations, sample-to-sample fluctu- 
ations can then be compared. Numerical studies have 
been carried out using different thermal-averaging meth- 
ods: the Metropolis Monte Carlo algorithm and exact 
enumeration. 

Due to the probabilistic nature of the Metropolis 
Monte Carlo algorithm, it is important to ensure that 
the system has reached equilibrium, and that the Monte 
Carlo average indeed gives the thermal average. Due to 
the long-range interaction term, we expect that a large 
system will need an extremely long time to reach equi- 
librium. At higher temperature, stronger disorder (i.e. 
larger W) and larger \/j,\, this difficulty may be avoided 
since the contribution of the Coulomb interaction is then 
small. It turns out that at low temperature fluctuations 
due to lack of equilibration are rather large up to the 
biggest(6x6) system we have reached. In order to study 
the fluctuations due to interactions in more detail, we 
turn to calculate the thermal average using exact enu- 
meration. In doing so, what occurs in lower temperature 
can be investigated more explicitly. 

Exact enumeration means that the thermal averages 
are achieved by finding explicitly all possible state config- 
urations and calculating their Boltzmann factors. From 
the comparison of results using two different methods, it 
is clear that those small-amplitude fluctuations obtained 
in using the Monte Carlo algorithm are due to the non- 
equilibration, instead of the interaction. 

With reasonable computing time, we obtain results for 
systems of sizes up to 5x4. Here, only the results for 
5x4 systems are presented. Fig. || clearly exhibits sig- 
nificant conductance fluctuations ranging from to e 2 /h 
as the chemical potential is varied. The amplitudes of 
the fluctuations, however, depend on the value of the 
energy scale Eq. We can choose different sequences of 
random numbers in the program to change the random 
site-energy at each site on the lattice. For the same dis- 
order strength (i.e. W), the sample-to-sample fluctua- 
tions due to different realisations of impurities can then 
be studied. In the two cases shown in Fig. there are 
large amplitude fluctuations and the details are sample 
dependent. 
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FIG. 2. Conductance versus chemical potential on a 5x4 
lattice at T=0.03, W=0.2, E =5, and h=0.3. The saddle 
point is chosen to be at (x,y)=(3,2). Two lines are for differ- 
ent realisations of impurity configuration. 

To avoid the confusions due to the arbitrary Eq, only 
plots of V s -p versus /i, instead of G versus /x, are shown in 
Fig. where h is the only free parameter to obtain V s — P . 
We then compare results at different disorder strengths. 
Fluctuations due to the Coulomb interaction should be 
suppressed as the strength of the disorder is increased. 
This is more or less consistent with what is demonstrated 
in Fig. ^. Finite temperature is another source in our 
model which can smear out the amplitudes of the fluctu- 
ations. This is because, with a fixed chemical potential, 
different thermal energies give different averaged occupa- 
tion numbers on the same lattice site. At higher temper- 
ature, each electron gains more thermal energy on aver- 
age and thus the effect of Coulomb interaction is again 
suppressed. The temperature dependence of the solid 
and dashed lines shown in Fig. ^ with the same disorder 
strength clearly shows this behaviour. 
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FIG. 3. Fluctuations of the saddle point potential, V B -p, 
at (x, y) — (3,2) versus chemical potential (or gate voltage) on 
5x4 systems, h =0.3 for all results here. 

Finally, we discuss the results on a system which con- 
tains randomly distributed localised states. This seems 
to be a more realistic representation of real devices. We 
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use the exact enumeration for the thermal averaging. The 
positions of the localised states are chosen independently 
with a uniform distribution over the system. For the 
same number of localised states as that in the lattice 
system, we conclude that finite-size effects are stronger 
than on the lattice. Although not shown here, results for 
different disorder strengths show no significant fluctua- 
tions. Comparing with the results for the lattice system, 
fluctuations are much smaller. In order to obtain promi- 
nent interaction-driven fluctuations off lattice, it seems 
possible that we need to study very large random-site 
systems. 

We now turn to describe the comparison between the 
numerical results from our model on a lattice and the 
experiments. First, let us focus on the conductance as 
a function of the gate voltage at some fixed magnetic 
field. The conductance fluctuations produced from our 
model are sample dependent, and varying the saddle- 
point potential with chemical potential can produce con- 
ductance fluctuations between G=0 and e 2 jh. Both of 
these features ace-consistent with what has been observed 
in experimentsoE3E3. However, the fluctuation patterns 
are rather different for our model and the experiments. 
It is observed in experiments that the conductance fluc- 
tuates in sharp peaks and-dips in the transition region 
between two plateauaij'EjO. By contrast, at least up to 
the system sizes for which we have done the computa- 
tion, it is observed that, from Fig. || and ||, sharp spikes 
of the fluctuations are obviously absent in our results. 
Instead, steps with rounded or flat tops or bottoms oc- 
cur as we vary /x. This feature can be understood in 
terms of the total electrostatic potential energy changing 
slowly or even remaining constant for some finite interval 
of the chemical potential. It is also true that the systems 
we have studied show a relatively small number of con- 
ductance fluctuations as /x varies, compared with experi- 
ments. However, we should expect more fluctuations in a 
bigger system. Our model seems to show this behaviour 
as the system size is increased. It is therefore not-appro- 
priate to analyse the distribution of conductancai-3 from 
our results which only contain small numbers of indepen- 
dent conductance values. 

Another important comparison is the temperature de- 
pendence of the conductance fluctuations. Finite tem- 
perature reduces the amplitudes of the fluctuations be- 
cause the fluctuation of the occupation number in each 
localised state is smoothed out by the finite thermal en- 
ergy. In our model, it can be observed that the width of 
each fluctuation 'steps' does not vary with temperature 
(see Fig. [|). In experiments, dataa demonstrate indeed 
that the amplitudes of peaks and dips are enhanced as 
the temperature is lowered. No obvious shrinking of their 
widths can be seen there. 

We now come to the behaviour of fluctuation peaks 
and dips at different magnetic fields. Although there is 
no explicit magnetic-field dependence in our model, we 
argue-here, by noting that straight lines on the V g -B 
planeca are obtained provided the system preserves the 



filling factor along each line, the occurrence of straight 
lines is consistent with the spirit of our model. 

As mentioned in the introduction section, the slopes of 
two sets of lines are parallel to those associated with the 
centres of the neighbouring plateaus. More explicitly, 
it is actually observed that each of these straight lines 
follows the equation V g = CB + D with C and D being 
constants. This behaviour can be, understood in terms of 
the physics of the Si MOSFETlJ, as has been discussed 
in Ref. 23. With the filling factor given by the relation 
v = (V g — V)esoh/de 2 B as a perpendicular magnetic field 
B is applied, we have 



Here V is some constant and represents a threshold volt- 
age, and d is a distance of the order of the thickness of 
the S1O2 layer. To describe the ith straight line on the 
Vg-B plot, another constant Vi, for example, is needed. 
For different parallel lines on the V g -B plane, we have dif- 
ferent Vi's. Along each line, the filling factor and V, are 
constants. This means that, as V g and B are both var- 
ied along the line, the fluctuation peaks and dips evolve 
in such a way that the filling factor of the system is un- 
changed for a given fluctuation. More explicitly, the value 
of v is observed to be either i or i + 1 for the (i + l)th 
transition region, depending on which plateau region the 
lines belong to. 

In our model, conductance fluctuations are associated 
with the occupancy of localised states in the 2DES. For 
the occupancy to fluctuate, the state must have energy 
near the chemical potential. Hence, along each straight 
line, which connects peaks or dips for different (V g ,B), 
the localised states are at the chemical potential. The 
fact that v along the straight line is the same as i, for ex- 
ample, associated with the ith plateau centre means that 
the localised states in both cases must have the same to- 
tal (kinetic plus electrostatic) energy as the states in the 
ith Landau level. This is indeed possible if the localised 
states belong to the Landau levels of the the 2DES. Fol- 
lowing this argument, as the chemical potential crosses 
the centre of a disordered-broadened Landau level in the 
transition region, localised states at two tails of Landau 
level then give the two slopes of straight lines correspond- 
ing to two different energies of the states at two plateau 
centres. 



IV. CONCLUSIONS 

In this paper, we have constructed a simple semiclas- 
sical model which produces conductance fluctuations in 
strong magnetic fields due to the Coulomb interaction be- 
tween electrons. The inclusion of the long-range Coulomb 
interaction causes the occupation number in each lo- 
calised state in the bulk to fluctuate as the chemical po- 
tential is varied. By taking into account the influence of 
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the electrostatic potential between localised states and 
the saddle point, the energy of the saddle point in the 
potential seen by mobile electrons then also varies with 
the chemical potential. Through this saddle point, the 
conductance due to the transmission of edge states from 
one contact to the other thus fluctuates with the chemi- 
cal potential. We study the model at finite temperature 
by numerical simulation using the Monte Carlo methods 
and exact enumeration. At low temperature, the Monte 
Carlo results suffer slow equilibration. Strong fluctua- 
tions due to the non-equilibration of the system in this 
case make it difficult to extract interaction-driven fluctu- 
ations. By contrast, results obtained using exact enumer- 
ation clearly exhibit significant fluctuations as a function 
of the chemical potential. 

In comparing these results with the experiments, our 
model shows qualitatively consistent behaviours with the 
experiments as the gate voltage and the magnetic field 
are both varied. There are, however, some different fea- 
tures existed between our results and experiments. In 
particular, although our simulations indeed exhibit fluc- 
tuations depending on realisations of disorder they give 
fluctuations which are like steps instead of the sharp 
peaks and dips observed in experiments. These discrep- 
ancies could arise because we have neglected interactions 
between bulk extended electrons and many specific de- 
tails in the Si MOSFET. 
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